Magnetic circular dichroism from the impurity band in III-V diluted magnetic 

semiconductors 
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The magnetic circular dichroism of III-V diluted magnetic semiconductors, calculated within a 
theoretical framework suitable for highly disordered materials, is shown to be dominated by optical 
transitions between the bulk bands and an impurity band formed from magnetic dopant states. 
The theoretical framework incorporates real-space Green's functions to properly incorporate spatial 
correlations in the disordered conduction band and valence band electronic structure, and includes 
extended and localized electronic states on an equal basis. Our findings reconcile unusual trends in 
the experimental magnetic circular dichroism in III-V DMSs with the antiferromagnetic p-d exchange 
interaction between a magnetic dopant spin and its host. 
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In III-V diluted magnetic semiconductors (DMSs) such 
as Gai_ x Mn x As, formed by doping nonmagnetic host 
semiconductors such as GaAs with magnetic acceptors 
such as Mn, the magnetic properties are highly corre- 
lated with the electrical and optical properties [U, H, H, H| • 
The effect of the added spin-polarized holes from the 
magnetic dopants on the relative optical absorption of 
right (<r + ) and left (c - ) circularly polarized light (mag- 
netic circular dichroism, or MCD) has been extensively 
explored to probe the interactions between the d states 
of the local magnetic moments and the p states of the 
host valence band @, [E 0, IE IE EE EH • Yet fundamen- 
tal puzzles remain: the light polarization most absorbed 
in Gai_ x Mn x As, which determines the sign of the MCD 
signal, is different from that of Zn!_ x Mn x Se. This sug- 
gests that the p-d exchange interaction in Gai_ x Mn x As 
is ferromagnetic [l0| . however in the dilute limit the in- 
teraction between a Mn spin and the GaAs valence band 
is known to be antiferromagnetic [13] • Refs. @, EH argue 
the unexpected sign of the MCD signal can be reconciled 
with an antiferromagnetic p-d interaction by considering 
a large shift of the Fermi level in the valence band due 
to doping (Moss-Burstein shift). However, the required 
Moss-Burstein shift is large (~ 100 meV) and the MCD 
spectrum should have a pronounced doping dependence 
that includes changing sign at low doping [J , whereas the 
observed doping dependence of the dominant features of 
the MCD spectrum is weak [6] and the unexpected sign 
of the MCD signal is present in low-doped, paramagnetic 
Gai_ x Mn x As (even for x ~ 0.005) 10]. Furthermore, 
none of the above treatments adequately treat the role 
of disorder in the optical transitions, whereas measure- 
ments indicate the carrier mean free path to be less than 
1 nm 14], comparable to the Fermi wavelength. Exper- 
imental evidence that the Fermi level in Gai-^Mn^As 
lies in the impurity band rather than in the valence band 

H EE E3, EE EE HE EH also requires reconsideration 



of the Moss-Burstein shift. 

Here we find that a proper consideration of the strong 
spatially-localized perturbation of the electronic struc- 
ture in the valence band due to the Mn dopant, including 
both the bound state of the acceptor and the perturba- 
tions of the continuum states near the dopant, is essential 
to accurately calculate the magnetic circular dichroism. 
A key element of the successful calculation of these prop- 
erties is our approach to impurity averaging the optical 
absorption in magnetic semiconductors. Traditional ap- 
proaches impurity average separately over the conduction 
band electronic structure and the valence band electronic 
structure, and then calculate the optical absorption for 
transitions between the two new effective bands as one 
would calculate optical absorption in a clean semiconduc- 
tor. Instead we calculate the difference in optical absorp- 
tion between the clean host and the host by calculating 
optical matrix elements between real-space Green's func- 
tions for a single dopant before impurity averaging. The 
effect of the short mean free path is included naturally 
by restricting the real-space sum to a small cluster of ap- 
proximately that diameter. Thus optical transitions that 
do not conserve crystal momentum are included with- 
out artificially relaxing momentum conservation between 
impurity-averaged bands 0] . 

In addition to obtaining the correct sign of the MCD 
signal for low doping as well as high doping, our calcu- 
lations identify the dominant transitions contributing to 
the MCD to be transitions between the bulk bands and 
the acceptor states bound to the Mn (that will form the 
impurity band). These transitions would be forbidden if 
impurity averaging and momentum conservation in opti- 
cal absorption calculations were improperly imposed, as 
both initial and final states would have different crys- 
tal momentum. We further find that the amplitude of 
the MCD tracks the magnetization of the material, as 
seen experimentally in the temperature-dependence of 
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FIG. 1: (color online) Absorption coefficient (ao) for GaAs 
as function of photon energy (fkv) with the Fermi level at 0.1 
eV (solid line), eV (dotted line), and —0.1 eV (dashed line) 
relative to the valence-band maximum. The step increase near 
1.9 eV shows the onset of the split-off band. The inset shows 
an example of Mo(E, E + Twj) as a function of E with fiu) — 2 
eV. The three peaks from highest to lowest correspond to the 
cases where the initial states are in the heavy-hole (HH), the 
light-hole (LH), and the split-off (SO) valence bands. 



the MCD [6(. Furthermore, the absorption onset is not 
sharp, and the MCD signal persists even for photon ener- 
gies below the band gap of the host semiconductor, also 
as seen experimentally 0| ■ 

We begin by describing the real-space tight-binding 
Green's function framework that permits the optical ab- 
sorption of a single dopant to be calculated without im- 
purity averaging. The real-space Green's functions to 
be used have already been used in the calculation of 
the electronic structure of Mn dopants and have shown 
excellent agreement with experimental measurements of 
the local electronic structure near individual Mn dopants 
H [H Hil and pairs of Mn [H [llembedded in GaAs. 

Following the formulation in Ref. |26| the absorption co- 
efficients for circularly polarized light (er^) can be written 
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dEf(E) [l-f(E')]M ± (E,E') ,(1) 



where E is the initial state energy, E' = E + hcu the final 
state energy, hio the photon energy, e the electron charge, 
n the index of refraction, Eq the vacuum permittivity, m 
the bare electron mass, c the speed of light in vacuum, 
f(E) the Fermi-Dirac function, and M ± (E,E') the opti- 
cal transition strength per unit energy per unit volume. 
In the electric dipole approximation 



M ± (E,E') = -L t r \piA(E')p ± A(E) 
mv L 



(2) 



where V is the system volume, p± = p x ± ip y are the 
momentum operators (the light is propagating along the 



z direction), and the spectral function, A(E), is a com- 
bination of retarded and advanced Green's functions, 



A{E) = - G R (E) 

Z7T 



G A (E) 



(3) 



The trace in Eq.([2]) is taken over a set of Lowdin orbitals, 
</>a,;( r — R-j,a), where j labels the primitive unit cells, 
a labels the atomic sites within a unit cell, Rj, a is the 
position vector of an atomic site, and I labels the atomic 
orbitals at each site. To simplify our calculations, we 
assume the momentum matrix elements are nonzero only 
between two Lowdin orbitals located at the same site and 
are independent of the type of the atom at that site. In 
the 16-band sp 3 tight-binding model we are left with only 
one type of momentum matrix element, 



(r-R^l^l^^r-Rj-a)) = iP . (4) 



With the above simplifications P can be linked directly to 
the momentum matrix element between the conduction 
and valence band Brillouin-zone-center states, known 
from bulk\i-p theory. Therefore, no additional parameter 
in our tight-binding framework needs to be introduced. 
The situation becomes more complex (and less empiri- 
cally constrained) if momentum matrix elements between 
neighboring atomic sites are permitted to be non-zero. 
Using the tight-binding parameters in Ref. |27J, the mo- 
mentum matrix element between the zone-center states is 
(T 4v<x \p x \T lc ) = 0.914 iP. Thus we find P 2 /m = 17.3 eV 
by setting the zone-center momentum matrix element to 
be the same as used in Ref. |28|. We calculate the tight- 
binding Green's functions with an energy linewidth of 
10 meV as in Ref. [22J. Lastly, we use a constant index of 
refraction, n = 3.878 for bulk GaAs at 2 eV ^ for the 
energy range (1.2-2.5 eV) shown in this Letter. 

We first apply this method to calculate the absorption 
coefficient for bulk GaAs. In this case Eq. ([2]) can be 
evaluated exactly in momentum space and the results, 
a ± = a and M ± — M , are shown in Fig. [TJ A good 
agreement with the experimental data (3pj as a function 
of the absorbed photon energy is found, although the 
calculated overall magnitude is larger by ~ 50 %, prob- 
ably due to the slightly too large conduction band mass 
typical for tight-binding models of III-V semiconductors. 
The energy dependence near the absorption edge has a 
square-root dependence with an Lorentzian tail because 
the quasiparticles have a finite lifetime and excitons are 
not included in our calculations. At higher energy the 
absorption's energy dependence is approximately linear. 
The shoulder near 1.9 eV is the onset for the split-off 
band (0.36 eV below the valence-band maximum). The 
absorption curves (solid and dotted lines) are approxi- 
mately the same as long as the Fermi level is in the gap. 
If the Fermi level lies in the valence band, the absorp- 
tion edge moves to higher energy and the shape is sig- 
nificantly altered. This shape change is not observed for 
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FIG. 2: (color online) Absorption coefficients, Aq + (dashed 
line), Aa~ (dash-doted line), and a~ — a + (solid line), as 
functions of photon energy, hui, with the Fermi level at (a) 
0.1 eV, (b) 0.05 eV, (c) eV, and (d) -0.1 eV relative to the 
valence-band maximum. 



Gai-a-Mn^As, and thus adds support to interpretations 
placing the Fermi level in the impurity band rather than 
in the valence band. 

In the presence of Mn dopants there are both local- 
ized states and extended states. Eq. @ is now evaluated 
in real space by summing over the lattice sites within 
in a finite cluster. We assume that the total MCD is 
the sum of MCD from these disjoint clusters contain- 
ing single Mn dopants. In doing so we have neglected 
Mn-Mn interactions. This approximation is sufficient 
to describe structures in the optical absorption that are 
broader in energy than the energy splittings due to Mn- 
Mn interactions (~ 0.1 eV for a Mn-Mn separation ~ 1 
nm [22j]). Calculations based on disjoint clusters con- 
taining single Mn dopants automatically produce MCD 
features that are directly proportional to the material's 
magnetization, consistent with the behavior of the prin- 
cipal experimental MCD features 0]. The cluster size is 
sufficiently large to account for transitions involving the 
localized Mn acceptor states or spin-polarized scattering 
resonances. The transitions between extended states ex- 
hibit finite-size effects in these cluster sums, but those 
do not contribute to MCD, and can be more accurately 
captured by the momentum-space calculation of Fig. [1] 
Therefore, only the difference AM ± = M ± — M is eval- 
uated using finite clusters, where M refers to the results 
without Mn. 

In our calculations, the light beam is parallel to one of 
the crystal axes and is also parallel to the Mn mag netic 
moment. The model of Mn is described in Ref. [22|. The 
acceptor states are bound mainly by the p-d exchange in- 
teraction, which is described by a spin-dependent poten- 
tial (V p d) present at the four lst-nearest-neighbor sites. 
V p d = 3.634 eV 24| is set to obtain the experimental 



binding energy (113 meV). An on-site potential (V n ) ac- 
counts for the direct Coulomb contribution to the binding 
energy, and is chosen to be 1 eV. Thus, 
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where (c) is the creation (annihilation) operator for 
electrons, n labels atomic sites (Mn is at n = 0), I labels 
atomic orbitals, and s labels spins. The spin quantization 
axis is parallel to the light propagation direction. The 
spin orientation of the Mn core 3d electrons is antiparallel 
to the light propagation direction, whereas the spin of 
the acceptor states is parallel to the Mn core spin. Thus 
Vpd > represents an antiferromagnetic p-d exchange 
interaction from the perspective of holes. The acceptor 
states are almost fully spin polarized and split into three 
energy levels due to the spin-orbit interaction[22]. As a 
result, the three levels have quite distinct orbital-angular- 
momentum character, and the top and bottom levels are 
coupled to opposite polarizations of circularly-polarized 
light. The upper level (farthest from the valence band 
edge) has orbital angular momentum parallel to spin. 

Our MCD results were obtained with a quasi-spherical 
cluster enclosing 99 atoms, with one Mn atom in the cen- 
ter. Up to the 8th nearest neighbors are included, which 
are 8 A away from the Mn. ~ 70% of the acceptor state 
is included in this cluster, which corresponds to approx- 
imately 2% Mn concentration. The results for different 
hole carrier concentrations are shown in Fig. [2] and we 
can see that the relative size of or is always larger than 
q + , no matter the doping. Note that we have ignored 
the orbital mixing in the impurity band by Mn-Mn in- 
teraction, which would reduce the MCD, and the energy 
broadening is taken into account only through a 10 meV 
linewidth in the single-particle Green's function. There- 
fore, the overall calculated MCD magnitude we obtain is 
larger than experimentally measured. 

To understand the nature of the doping-independcnt 



"positive" (a 



> a 



') MCD signal, the contributions 



from different optical transitions at the typical photon 
energy (2 eV) are shown in Fig. [3] as a function of the 
initial-state energy. Two peaks for each circularly polar- 
ized light are present, corresponding to the processes as- 
sociated with the transition of an electron to or from the 
impurity levels. In our configuration, the upper impurity 
level has predominately angular momentum projection 
— 1 and the bottom level has projection +1 because the 
spin-orbit interaction favors the orbital angular momen- 
tum parallel to spin (in the valence electron convention) . 
The processes between the valence and conduction band 
are indicated by the short arrow in Fig. [3j The sharp- 
ness of the peaks of M in the energy space comes from 
the small broadening factor (10 meV) that we used for 
the impurity level. From Fig. [3] it is apparent that both 
cr + processes are suppressed by the Pauli exclusion prin- 
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FIG. 3: (color online) The optical transition strength, 
AM+(E,E + ?uj) (dashed line), and AM~ (E, E + tkj) (solid 
line), as functions of E with huj = 2 eV. The four peaks cor- 
respond to the four distinct processes shown by the arrows in 
the energy diagram at the center. States in the shaded region 
are occupied if the Fermi level lies in the impurity band. 
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FIG. 4: (color online) The absorption-coefficient differences, 
Aa + (dashed line) , Aa - (dash-doted line), and a - — a + (solid 
line), with hui — 2 eV as a function of the Fermi level relative 
to the valence-band maximum. 



ciple if the Fermi level lies in the middle of the impurity 
band. This is further illustrated in Fig. 2] We see that 
the dichroism is the strongest when the Fermi level is 
between the upper and the lower impurity levels. The 
suppression of the positive MCD signal when the Fermi 
level is above the impurity band is consistent with recent 
observations in donor-compensated systems [ill ]. 

In conclusion, a theoretical approach to calculating the 
optical properties of DMSs using real-space tight-binding 
Green's functions naturally explains several of the puz- 
zling experimental findings that have been insufficiently 
explained using impurity-averaged conduction and va- 
lence states. The MCD measurements agree with anti- 
ferromagnetic p-d exchange interaction expected from 
single-Mn measurements, and no Moss-Burstein shift is 
required to explain the sign of the MCD measurements. 



There is no sharp absorption onset due to transitions 
from the valence band to the impurity band and the 
dominant MCD features are proportional to the material 
magnetization, as previously reported from experimental 
observations. 

This work was supported in part by an ONR MURI 
and an ARO MURI. 
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